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1 INTRODUCTION 



it has been shown that 
play a significant role in 



Through recent observations 
active galactic nuclei (AGN) 
the evolution of galaxies. The observed correlation be- 
tween the mass of the central black hole and the velocity 
dispersion of the bulge of the host galaxy suggests a 
strong connecti on between galaxy evo l ution and black hole 



ABSTRACT 

Using a fully cosmological hydrodynamic simulation that self-consistently incorporates 
the growth and feedback of supermassive black holes and the physics of galaxy for- 
mation, we examine the effects of environmental factors (e.g., local gas density, black 
hole feedback) on the halo occupation distribution of low luminosity active galactic 
nuclei (AGN). We decompose the mean occupation function into central and satellite 
contribution and compute the conditional luminosity functions (CLF). The CLF of the 
central AGN follows a log-normal distribution with the mean increasing and scatter 
decreasing with increasing redshifts. We analyze the light curves of individual AGN 
and show that the peak luminosity of the AGN has a tighter correlation with halo 
mass compared to instantaneous luminosity. We also compute the CLF of satellite 
AGN at a given central AGN luminosity. We do not see any significant correlation 
between the number of satellites with the luminosity of the central AGN at a fixed 
halo mass. We also show that for a sample of AGN with luminosity above 10 42 ergs/s 
the mean occupation function can be modeled as a softened step function for central 
AGN and a power law for the satellite population. The radial distribution of AGN 
inside halos follows a power law at all redshifts with a mean index of —2.33 ± 0.08. 
Incorporating the environmental dependence of supermassive black hole accretion and 
feedback, our formalism provides a theoretical tool for interpreting current and future 
measurements of AGN clustering. 

the black hole mass fu nction (e.g., IShankar et al. 2004: 
I Graham 



hole mass fu nction (e.j_ 
Driverl 120071 ) measured at the current epoch 



activ ity 

l200ll ; iTremaine et al 



Gebhardt et al 
2002T T 



|2000|; iMerritt fc Ferraresel 



can provide constraints on models of black hole growth. 
Measurement of AGN clustering provides a unique way 
to study the physical characteristics of AGN through 
the connecti o n wit h their hos t dark matter halos (e.g. , 
Croom et al] |2004|; 



Shen et al . 2009; R oss et all [2009) 



Gilli et all 120051; iMvers et al.l 120061 ; 



Graham et alj 



20 111 ). 



Sev- 
eral theoretical models relating AGN activity or black 
hole growth to g alaxy evolution have been proposed 



gro 

(e.g.. ISoitanl 119821; ISilk 



Kauffmann fc HaehneltJ 



2006 



2006 



'z Rees 1998; Salucci ct al. 
l2000l; IWvithe fc Loebl 



Marconi et~all |2004|; JCattaneo et all 120061; ICroton et al 



DTMatteo et ai1l2005l ; lHopkins et al.ll2006l ; lLapi et al 
IShankar et al. r i2004l ). Differentiating between these 
theoretical models requires several observational quantities. 

Measurements of the AG N luminosity function (e.g., 

iBovle et al1l200d ; iFan et al.ll200ll ) and the number density 
of black hole hosts in the present universe can provide an 
estimate of the duty cycle of black holes. Alternatively, 



Theoretically clustering properties of AGN have been 
mostly studied with semi-analytic models using the halo 
model or the black h o le continuity equat i on approach (e.g . , 
IShankar etabl l2010al ; iBonoli et all 120091 ; iLidz et all |2006| ). 
Although semi-analytic models have provided the formal- 
ism for interpreting quasar clustering, there are certain is- 
sues that these models cannot approach. In this formalism 
physical parameters (e.g., "quasar duty cycle") are treated 
as free paramete rs that are constrained from clustering mea- 
surements (e.g.. [Martini fc W einberg 2 00 1] ). There also ex- 
ist degeneracies betw een parameters in the models (e.g., 
IShankar etabl l2010bl ). Moreover in most of these models 
black holes are accreting at a fixed fraction of the Eddington 
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Boxsize N p 


7UDM 


'^gas 


e 


{h' 1 Mpc) 






(h- 1 kpc) 


33.75 2 x 486 3 


2.75 x 10 7 


4.24 X 10 6 


2.73 



Table 1. The numerical parameters in the simulation. N p , ihdm, 
m g as and e are defined as the total number of particles, mass of 
the dark matter particles, mass of the gas particles, and comoving 
gravitational softening length respectively. 



rate, making accretion rate dependent on black hole mass 
only. In reality, parameters like accretion efficiency and the 
duty cycle are not fixed and depend strongly on environ- 
ment (e.g., local gas density, mergers, feedback from AGN). 
Some recent semi-analytic studies started to consider more 
general cases of black hole accretion to investigate the cos- 
mological co-evolution of black holes and the evolu t ion of 
AGN luminosity functions (e.g.. [Marulli et ail |200S| , |2009| ; 
iMalbon et al.ll2007l ; lMenci et al.ll2008l ). 

Hydrodynamic simulations of galaxy formation with 
black hole growth can capture the environmental depen- 
dence of accretion. The growth of black holes is tied with 
the full dynamics of dark matter. Also with these sim- 
ulations we can model the interplay between AGN and 
galaxy form ation in the form of fee dback and self-regulated 
growth fe.g.. lWvithe fc Loebl l2003h . Thus these simulations 
provide an excellent platform to study the co-evolution of 
AGN with large scale structures over cosmological time 
scales. Studies of AGN clusterin g using cosmological sim- 
ulations have been ca r ried out by[ Thackcr et al. (2009J) and 
iDegraf et all (|2011al ). iDegraf et all (|2011al ) calculates the 
correlation function of black holes using a smoothed particle 
hydrodynamics (SPH) cosmological simulation that incorpo- 
rates galaxy formation physics and self-consist ent black hole 
growth and feedback (|Di Matteo et all l2008h . They show 
that the black hole correlation function consists of two dis- 
tinct components: contributions from intra-halo and inter- 
halo black hole pairs, i.e., the one-halo and two-halo terms. 
At small scales the one-halo term still follows a power law, 
which is different from that of dark matter. This boost in 
small scale power is due to galaxies hosting multiple black 
holes. 

Despite their advantages, studies using hydrodynam- 
ical simulations are limited due to their computational 
costs. For examp l e sim ulations are limited to small boxes 
|Di Matteo et all l200ct ) restricting the analysis to small 
scales and low luminosity AGN. The computational cost of 
running these simulations to z = limits the possibility 
to compare the outputs with observations in the local uni- 
verse. Also since every modification of the recipies adopted 
to describe the 'sub-grid' physics requires rerun of the simu- 
lations, it is not feasible to do systematic parameter studies. 
Semi-analytic studies are advantageous over hydro simula- 
tions in this context. W e now employ the SPH simulation of 
|Pi Matteo et all (|2008l ) to investigate the relation between 
AGN and their host halos within the halo occupation dis- 
tribution (HOD) framework, a useful analytic formalism for 
modeling and interpreting AGN clustering. In our approach, 
the expensive hydrodynamic simulation is used to obtain in- 
sights on a general analytic technique for studying AGN 



clustering. The HOD (e. g., iMa fc Fry! l200d ; ISeliaJj l200d : 
iBerlind fc Weinberg! I2OO2T) provides an analytic formalism 
for understanding clustering properties of galaxies or, in this 
case, AGN. It is characterized by the probability P(N\M) 
that a halo of mass M contains N AGN of a given type, 
together with their spatial and velocity distribution inside 
halos. As long as the HOD at a fixed halo mass is statisti- 
cally independent of t he large scale environments of halos 
(e.g., iBond et al.lll99ll ). it provides a complete description 
of the relation between AGN and halos, allowing the calcu- 
lation of any clustering statistics (e.g., two point correlation 
function, three point correlation function, void probabil- 
ity distribution, pairwise velocity distribution) at all scales 
(small, intermediate, large) for a given cosmological model. 
This implies that if we can constrain the HOD empirically, 
we will have knowledge of everything that the measured 
clustering properties have to tell us about AGN formation 
and evolution models. The HOD allows to distinguish be- 
tween the background cosmology and AGN evolution mod- 
els. The cosmological parameters are encoded in the distri- 
bution of halos, where as the bias between mass and AGN is 
fully described in the probability distribution P(N\M ) (see 
IBerlind fe Weinberg! I2OO2I IZheng fc Weinberg! 120071 for the 
strength of HOD). The HOD formalis m has been widely u sed 
in interpreting galaxy clustering (see lZehavi et al]|2005l ). 

In this paper we perform a theoretical study on the rela- 
tion between AGN and dark matter halos using cosmological 
simulations. We assess the validity of several simplifying as- 
sumptions in the HOD modeling of AGN and its predictive 
power and limitations for modeling AGN clustering data. 
Our paper is organized as follows: In §2 we describe our 
simulations and give a brief description of our methodology. 
In §3 we present the conditional luminosity function of AGN, 
the mean occupation distribution, and the radial profiles of 
AGN. Finally in §4 we summarize our main conclusions and 
provide future directions to this work. 



2 SIMULATION 
2.1 Numerical Code 

The numerical code uses a standard ACDM cosmological 
model with cos mological paramete rs from the first year 
WMAP results jSpergel et alj|2003h . The simulation uses 
an extended version of the paralle l cosmological Tree Par- 
ticle Mesh-SPH code GADGET2 dSpringe]||2005l). Gas dy - 
namics is modeled with Lagrangian SPH (|Monaghan| [l992'): 
radiative cooling and heat ing processes are computed from 
the prescription given by iKatz et all (|l99d ). The relevant 
physics of star formation and the associated supernova feed- 
back have been approximated based on a sub-resolution 
multi phase model for th e inte rstellar medium developed 
by ISpringel fc HernquistJ |2003l ). The size of the simula- 
tion box is 33.75/i -1 comoving Mpc with periodic bound- 
ary conditions. A detailed description of the implementa- 
tion of black hol e accretion and t he as sociated feedback 
model is given in |Pi Matteo et all (|200ct ) . Black holes are 
represented as collisionless "sink" particles that can grow in 
mass b y accreting gas or by mergers. The Bondi-Hoyle re- 
lation (|Bondilll952h is used to model the accretion rate of 
gas onto a black hole and capture the environmental depen- 
dence of black hole accretion. The accretion rate is given 
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Figure 1. Distribution of accretion rates scaled by the Eddington value (\ = M / Me~d~d) as a function of halo mass at redshift 5.0 (left 
panel), 3.0 (middle panel), and 1.0 (right panel). The black points show the distribution of the central AGN while the red open circles 
show the distribution of the satellite AGN. The mean rate depends weakly on halo mass at all redshifts. For a given halo mass, AGN at 
higher redshifts accrete at higher rates than at low redshifts. 



by Mbh = 47r[G 2 M| H p]/(Cs + v 2 ) 3/2 , where p and c s are 
density and speed of sound of the local gas, v is the velocity 
of the black hole with respect to the gas, and G is gravi- 
tational constant. Although the Bondi parameterization as- 
sumes spherical accretion, observations show that the Bondi 
relation adequately captures the physical state of the black 
hole and the Bondi sca ling holds at scale s much larger than 
the Bondi radius (e.g.. I Allen et al.ll2006l ). 

The bolometric luminosities of the AGN are LboI = 
»?A4bhc 2 where r\ — 0.1 is the canonical efficiency for thin 
disk accretion. It is assumed that a small fraction of the 
bolometric luminosity couples to the surrounding gas as 
feedback energy Ef such that Ef = e/LboI with the feed- 
back efficiency e/ taken to be 5%. Based on previous stud- 
ies of galaxy mergers, the value of the feedback efficiency 
parameter e/ is chosen to be 5% in the simulation, a nec- 
essary fraction to reproduce the n ormalization of the ob - 
served Mbh — cr relation at z — (|Di Matteo et all l2005h . 
This feedback energy is put directly into the gas smooth- 
ing k ernel at the position of the black hole |Pi Matteo et all 
2008). The feedback energy is assumed to be distributed 
isotropically for the sake of simplicit y; however the response 
of the gas can be ani sotropic (e.g.. IChatteriee et al.l 120081 ; 
|Pi Matteo et al]|2008l ). The formation mechanism for the 
seed black holes which evolve into the observed supermassive 
black holes today is not known. The simulation creates seed 
black holes in halos that cross a specified mass threshold. 
At a given redshift, halos are defined by a friends-of-friends 
group finder algorithm run on the fly. For any halo with 
mass M > 5 x 10 10 /i -1 Mq that does not contain a black 
hole, the densest gas particle is converted to a black hole of 
mass Mbh = 5 x 10 s /i -1 Mq. The choice of the seed mass in 
the present simulation was based on curren t models of seed 
black hole formation in the u niverse (e.g.. iBromm fc Loebl 
120031 ; IBromm fc Larsonl 1200 j ). The black hole then grows 
via the accretion presc ription given above an d by mergers 
with other black holes l|Di Matteo et al.|[2008l ). The simula- 
tion parameters are shown in Table 1. 

The detailed parameter studies of this simulation, 
and comparison with several observations have been pre- 



viously done by Siiacki et alj (12007 ) ; Pi Matteo et all 
j200St); IColberg fc Pi Matted (|2008t ); fcroft et all |2009l ); 
iPegraf et al l (|2010l ). The model has reproduced the ob- 
served Mbh — cr relat ion, total black hole mass density 
l|Pi Matteo et al1l2008l ). and the AGN luminosity functions 
and their evolution in optical, soft and hard X-ray band 
jPegraf et al.ll20ld) . The black hole mass density from the 
simulation matches well with the con straints from the in- 
tegrated X-ray lum inosity function (|Shankar et al.l I2004J : 
M arconi et al.l 12004 ) and the accretion rate density is con- 
sistent with the constraints of IHopkins et all (|2007al ). This 
model has been also used in galaxy merger simulations to 
investigate bl ack hole growth and their correlation with 
host galaxies |Pi Matteo et al.ll2005l ; iRobertson et al.ll2006l ; 
IHopkins et al.ll2007bl )~ 

The model has some intrinsic limitations. It does not 
treat the physics of the accretion disc in detail, which is 
not feasible in a simulation with cosmological volume. Nev- 
ertheless, the model is capable of representing some key 
aspects of black hole evolution in a cosmological context 
and in reproducing many observational results. Although 
the model has some limitations, it has been fairly success- 
ful in reproducing some of the key observational results. 
Several other teams ( Booth & Schaye 2009; Joha nsson et al.l 
|200S| ; iTevssier et all I l201ll ) have now implemented similar 
modeling for black hole accretion and feed back in hydro- 
dynamic simulations. iBooth fc Schavei (2009) have indepen- 
de ntly explored the parameter space of the fidu cial model 
of iPi Matteo et all (|2008l ). ITevssier et ail (|201ll ) have im- 
plemented this prescription in an adaptive mesh refinement 
code (different from the SPH formalism) and showed that 
AGN feedback can solve the overcooling problem in clus- 
ters. This large number of previous investigations make this 
particular subgrid model a good choice for studying co- 
evolution of AGN with their host dark matter halos. We 
emphasize that in this work we do not intend to explore the 
parameter space related to the subgrid modeling of black 
holes. Instead we study the relation between AGN and dark 
matter halos in one simulation box. The goal of our cur- 
rent study is to obtain a useful description between AGN 
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and dark matter to guide the interpretation of observations 
(e.g., AGN cluster ing) and to learn abou t AGN evolution 
using the model of |Pi Matteo et all (|2008l ). 



2.2 Sample Selection 

For our analysis we consider halos above a mass scale of 
1O 11 M0 and black holes with masses above 1O 6 M0. This 
choice of mass scales will minimize the effect of the seed mass 
in our results. We also restricted our analysis to black holes 
with Eddington scaled accretion rate \ = M/Medd > 1CP 5 
since we want to study systems with active accretion. This 
selection imposes a luminosity cut of 10 38 ergs/s in our sam- 
ple. One of the key goals of this paper is to separate the con- 
tributions to the occupation distribution fro m the central 
and sa tellite AGN (as done for galaxy HOD bv lZheng et alj 
(2005)). We consider all the black holes within R200 (defined 
as the radius within which the enclosed mean density is 200 
times the critical density) of a halo to be associated with 
the host halo. The most massive black hole within -R200 is 
assumed to be the central AGN and the rest are designated 
as satellites. At z = 1.0,3.0, and 5.0 about 79%, 93%, and 
91% of the central black holes lie within 20% of -R200 respec- 
tively. We expect the central black hole to be located at the 
center of the main galaxy in the halo since this is the place 
of highest gas density and thus maximal black hole growth. 
The fact that most 'central AGN' are found near the group 
center shows that our technique for identifying the central 
AGN is reasonably adequate. Due to the limited volume of 
our simulation we co uld only probe the faint end of the AGN 
luminosity function (|Degraf et al-l feoiO), For calculating the 
mean occupation we study AGN samples with bolometric 
luminosities above 10 42 ergs/s, corresponding to the obser- 
vational limit of X-ray selected AGN in deep surveys (e.g., 
lUeda et al.|[2003l ; iBovle et al.lll993l ; TCowie et al.ll2003l ). We 
also show the distribution of AGN with luminosities above 
10 40 ergs/s and 10 38 ergs/s for theoretical interests to see 
the luminosity dependence of HOD parameters. 

In a companion paper, iDegraf et all (l2011bh (hereafter 
Paper I), we use the same SPH simulation to characterize 
the HOD of black holes as a function of black hole mass. 
Our present study is fundamentally different from that in 
Paper I. We aim to study the HOD of AGN in terms of 
AGN luminosity. By definition the HOD is always stud- 
ied based on physical properties of the object. In the case 
of galaxy clustering, the galaxy HOD has been extensively 
studied i n the literature based on color, lumin osity, stellar 
mass (see lZehavi et al.ll2005l . IZheng et al.ll2007l ). Black hole 
mass and luminosity (accretion rates) are different param- 
eters and there does not exist an obvious one-to-one corre- 
spondence between mass and accretion rates (e.g. Fig. 1). 
Additional baryonic physics is needed to go from black hole 
mass to accretion rate (and thus luminosity). There is some 
correlation between mass and accretion rates, but the scat- 
ter is huge and the correlation depends strongly on redshift 
l|Colberg fc Pi Matted l200sl ; IChatteriee et alj|2008h . 

In practice, it is barely possible to have accurate mass 
measurement of a large sample of black holes beyond the 
local universe. Therefore, the study in Paper I is mostly of 
theoretical interest. By investigating the occupation prop- 
erties of black holes at the subhalo level, Paper I develops 
an analytic mechanism to populate dark matter halos with 



black holes in simulation. On the contrary, AGN luminos- 
ity is more readily measured observationally and there are 
existing and ongoing efforts in measuring AGN/quasar clus- 
tering as a function of luminosity. Our study in this paper 
will provide the necessary ingredient in interpreting these 
observations. Similar to the mass function of black holes 
and luminosity function of AGN, the black-hole-mass based 
HOD and AGN-luminosity based HOD are complementary 
tools to advance our understanding of black hole evolution. 

Since the goal of this paper is to provide the HOD 
framework to interpret observations we adopt some tech- 
nical differences in the definition of dark matter halos and 
the choice of our sample from Paper I. In Paper I we dif- 
ferentiate between central and satellite black holes at the 
subhalo level. Black holes residing in the central subgroup 
are called central black holes and black holes residing in 
satellite subgroups are called satellite black holes. In the 
nomenclature in Paper I there can be multiple central black 
holes. We emphasize that the model in Paper I will be useful 
for populating halos (N body simulations) with black holes 
using semi-analytic approaches and for studying the distri- 
bution of black holes in cosmological simulations. On the 
other hand our approach in this paper is to perform a the- 
oretical study on the relation between AGN (selected based 
on luminosity) and dark matter halos and examine the as- 
sumptions in the HOD for modeling AGN clustering data. 
Our identification of central (i.e. most massive black hole 
within -R200) and satellite (rest of the black holes within 
R200) AGN is equivalent to primary and secondary black 
holes in Paper I (see Fig. 2 of Paper I) and is similar to the 
standard terminology used for galaxy HOD. We note that 
the friends-of-friends halos were used in Paper I for study- 
ing black hole occupation. Here we choose -R200 to define 
the halo boundary, which is commonly adopted in defining 
galaxy clusters. Since the HOD describes bias at the level of 
systems near dynamical equilibrium, traditional virial meth- 
ods of mass estimation can take advantage of this and pro- 
vide a di rect measurement of the p robability distribution 
P(N\M) (|Berlind fc Weinberd[2002T ). We have verified that 
the two different halo definitions show a good agreement in 
assigning black holes to host halos within 10%. 



3 RESULTS 

In this section we present the distribution of accretion rates 
as a function of halo mass and black hole mass. From this 
we calculate the conditional luminosity functions and the 
mean occupation distribution of central and satellite AGN. 
Finally we show the radial distribution of the satellite AGN 
and provide the best-fit hod parameters. 



3.1 The Accretion Rates 

Figure 1 shows the distribution of the Eddington scaled ac- 
cretion rate (x = M/Medd), as a function of halo mass 
for three redshifts. The black points are central AGN and 
the red open circles represent satellites. The mean rate is 
roughly independent of halo mass, which is commonly as- 
sumed in semi-analytic studies. However, the mean value of 
X at a given halo mass varies with redshift. AGN at higher 
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Figure 2. The distribution of Eddington scaled accretion rates (x) for different black hole mass and redshift bins. From left to right, 
columns correspond to black hole mass bins, 6.0 < Log(MeH) < 6.5, 6.5 < Log(MBfr) < 7.0, and Log(MeH) > 7.0. The top, middle, 
and bottom panels correspond to redshifts 1.0, 3.0, and 5.0 respectively. The solid and dotted lines represent the central and satellite 
distributions. The mean of the distribution is lower at Z = 1.0 compared to redshifts 3.0 and 5.0 for all mass bins. The error bars are 
Poisson error bars. 



redshifts (left panel of Fig. 1 showing z = 5.0) tend to ac- 
crete at a higher value than their low redshift counterparts 
(right panel of Fig. 1 showing z = 1.0). Figure 2 shows 
the distribution of x as a function of black hole mass. It is 
defined as the conditional probability distribution of \ f° r 
a given black hole mass P(x\Mbu)- The solid and the dot- 
ted lines denote central and satellite AGN, respectively. The 
mean rates are again roughly independent of black hole mass 
but vary substantially with redshift. The central and satel- 
lite distributions tend to trace each other at all redshifts. 

In Fig. 3 we present the spatial distribution of the satel- 
lite AGN. The sizes of the circles are proportional to the 
black hole mass (left panel), gas mass (middle panel), and 
stellar mass (right panel) respectively. The gas and the stel- 
lar mass are computed within a spherical region of radius 
10 kpc surrounding the black hole. The green triangles rep- 
resent black holes with zero stellar mass within 10 kpc. The 
top, middle, and bottom panels are for z = TO, 3.0, and 
5.0 respectively. The AGN accretion rates do not show any 
dependence on radius or black hole mass. Also we do not 
see any correlation between black hole mass and its loca- 
tion within the halo. The trends of accretion rates with gas 
and stellar mass are more prominent. The low values of x 
at low redshift do not correlate with the mass distribution 
of the black holes (left panel of Fig. 3). The mass distribu- 
tion of black holes is not significantly different at different 



redshifts. At low redshifts the gas mass around black holes 
decreases and so do the accretion rates. This is due to the 
effect of AGN feedback. Outflow from AGN is responsible 
for expelling gas from the vicinity of the black hole and sup- 
pressing its own growth. AGN feedback is responsible for 
shutting down star-formation too. In Fig. 3 we also see the 
evidence of suppressed stellar mass in the immediate neigh- 
borhood of the AGN. We see evidence of how a local phe- 
nomenon (feedback) affects the global distribution of AGN, 
regardless of the host halo mass (see Fig. 1 for the host halo 
masses) . This picture is consistent with what has been seen 
in pr evious studies (e.g.. |Pi Matteo et a l. 2008; S iiacki et alj 
l2007h and semi-analytic models (e.g., IWvithe fc Loebll2003l ) 



3.2 The Conditional Luminosity Function 

The c onditional luminosity function (CLF) fe.g. jYang et alj 
2003), in this case for AGN, is defined as the luminosity 
distribution <&(L\M) of AGN that reside in halos of a given 
mass M. The global luminosity function is given by 



*(L) 



dn 

Jm 



${L\M)dM, 



(1) 



where $(£) is the AGN luminosity function and dn/dM 
is the halo mass function. The CLF constitutes the differ- 
ential form of the HOD. The CLF provides a tool to ex- 
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Figure 3. Spatial distribution of the satellite AGN. The size of each circle is proportional to the mass of the black hole (left panel), gas 
mass (middle panel), and stellar mass (right panel). The gas and the stellar mass are computed within a spherical region of 10 kpc around 
the black hole. The green triangles represent black holes with no stars within 10 kpc. The top, middle and bottom panels correspond to 
z = 1.0, z = 3.0, and z = 5.0 respectively. The distribution has been shown for all satellite AGN residing in halos with M^aio > lO n M0 
and black hole mass Mbh > 10 6 Mq. The accretion rates do not show any dependence on radius. The accretion rates do not show any 
pattern with black hole mass either, but there is a correlation with gas mass and stellar mass. In general the accretion rate is higher for 
a higher gas (stellar) mass. As we go to low redshifts the gas mass and stellar mass around black holes decreases due to the effect of 
AGN feedback. AGN feedback is responsible for expelling gas from the vicinity of the black hole and suppressing its own growth. 



amine the distribution of halo mass for a given luminos- 
ity and study luminosity dependent clustering. This for- 
malism has been widely used in modeling galaxy cluster- 
ing with data from 2dFGRS an d DEEP2 redshift surveys 
(e.g.. Ivan den Bosch et al"]|2003l ). The CLF is useful in gen- 
erating mock catalogs w hich can then b e used to test th e 
cosmological model (e.g., iMo et al ] |2004l : lYan et alJliool ). 
Figure 4 shows the CLF in bins of halo mass running from 
Log M = 11.25 ± 0.25 (left most panel) to Log M > 12.5 
(right most panel). The dashed lines show the distribution 
of bolometric luminosities for the central AGN while the 
dotted lines show the distributions for the satellite AGN. In 
both cases the distribution has been normalized by the to- 
tal number of halos in the given mass bin. The top, middle, 
and bottom panels show results for redshifts 1.0, 3.0, and 
5.0, respectively. 

The luminosity distributions reveal several features. 
The central AGN luminosity distribution closely traces a 
log- normal distribution (e.g., iMartini fc Weinberg! 1200 ll ). 
Fig. 5 captures the features we observe in the CLF. The 
top panel of Fig. 5 shows the evolution of the mean log lu- 
minosity of the central AGN as a function of halo mass. 



The bottom panel of Fig. 5 shows the la scatter in the log 
of the luminosity. The mean luminosity at a given halo mass 
evolves with redshift, with higher mean luminosity at higher 
redshifts. The scatter in the luminosity distribution of cen- 
tral AGN increases toward lower redshift, indicating that the 
luminosities for a given halo mass are more uniform at high 
redshift than at low redshift. As a result of this, at low red- 
shifts AGN samples based on luminosity will have a wider 
range of host halo masses and hence clustering will depend 
weakly on luminosity. However at high redshifts luminosity 
dependent clustering will be more prominent. Similar red- 
shift evolution of luminosity d ependent clusteri n g has been 
observed with SDSS quasars ijShen et al.l 120091 . 120071 ). We 
note that the clustering properties of AGN will be largely 
related to their host halo mass. However we find that ac- 
cretion rates and hence luminosity depends strongly on lo- 
cal properties particularly at low redshifts (e.g., feedback 
discussed in §3.1) which erases some of the dependences of 
bolometric luminosity on halo mass and hence the luminos- 
ity dependent clustering at low redshifts. This will not be 
the case in black hole mass selected sample since black hole 
mass is more tightly correlated with halo mass regardless 
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Figure 4. The conditional luminosity functions of the AGN in the simulation. The top, middle and bottom panels show the probability 
distributions for redshifts 1.0, 3.0, and 5.0, respectively. The halo mass bins range from 10 n M Q - lO n - 5 M (leftmost), to > 1O 12 - 5 M 
(rightmost). The dashed and dotted lines represent the probability distribution for the central and satellite AGN. The curves are 
normalized by the total number of halos in each mass bin. The central AGN traces a log-normal distribution as commonly assumed in 
semi-analytic studies. The error bars represent the Poisson error bars in each bin. 



of redshifts (|Colberg fc Pi Matted l2008h . The shape of the 
satellite luminosity function cannot be identified definitively 
due to lack of statistics. We note that the lower end of the 
satellite luminosity function is affected by the resolution of 
the simulation and the seed black hole mass, which is mani- 
fested in the cut-off of the luminosity function. The satellite 
distribution also shows some variation with halo mass. For a 
given halo mass the peak of the satellite distribution tends 
to be at a lower luminosity than the central AGN. 

To examine the effect of the central AGN on the num- 
ber distribution of satellites within a halo, we compare the 
conditional luminosity functions of satellite AGN for halos 
differing in central AGN luminosities. This is defined as the 
conditional distribution of satellite AGN luminosities for a 
fixed Mhaio and LbS- We used a large halo mass bin to in- 
crease the statistics of our sample. We note that there is a 
correlation between central AGN luminosity and halo mass. 
To eliminate the effect of mass-dependent central AGN lu- 
minosity, we divide halos according to the central AGN lu- 
minosity as follows. For each redshift, at each halo mass, we 
tag halos as 'high Lcen' ('low Lcen') if their central AGN 
luminosities are above (below) the mean central AGN lumi- 
nosity at that mass (Fig.5). Then the conditional luminosity 
functions of satellite AGN are computed for the 'high Lcen' 
and 'low Lcen' halos, respectively. The results are shown 
in Fig. 6. The solid lines show the distribution of satellite 



luminosities for 'high Lcen' sample and the dashed lines rep- 
resent the distribution of satellite luminosities for 'low Lcen' 
sample. 

The two groups of curves are for two halo mass bins. The 
open and the filled black circles connected by the dashed and 
solid lines respectively represent the halo mass bin 11.0 < 
Log(Mhalo) < 12.0. Similarly the blue open triangles and the 
open squares connected by the dashed and solid lines respec- 
tively represent the halo mass bin 12.0 < Log(Mhaio) < 13.0. 
The curves are normalized to the total number of halos host- 
ing the central AGN above and below the mean luminosity. 
We have shown the result for z = 3.0. Our results are similar 
for z = 1.0 and z = 5.0. The correlation between central and 
satellite AGN luminosity at a fixed halo mass will have im- 
portant implications on the small scale clustering strength. 
If there exists a strong large scale feedback effect from the 
central AGN it can possibly alter the local gas distribution 
around satellite AGN and suppress the growth of the satel- 
lite black holes. This in effect will decrease the number of 
luminous satellites in a halo with a higher central AGN lu- 
minosity and we would observe an anti-correlation between 
satellite number and central AGN luminosity. 

In Fig. 6 we do not see any correlation between the 
number of satellites a nd the luminosity of th e central AGN. 
It has been shown in lChatteriee et all (|2008l ) that m group 
scale halos (the most massive halos in our simulation) feed- 
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Figure 5. The mean (upper panel) and scatter (lower panel) of 
the central AGN conditional luminosity function as a function of 
halo mass, shown in Figure 4. The scatter is defined as the la 
scatter from the mean in LogLBol- The dot-dashed, solid, and 
dashed lines denote redshifts 1.0, 3.0, and 5.0, respectively. For 
a given halo mass the mean increases with redshifts. The growth 
of the black holes is possibly suppressed and hence the accretion 
rate (and the luminosity thereof) is lower at lower redshifts. The 
scatter in the distribution is also high at low redshifts since several 
factors (e.g., feedback from black holes) introduce spread in the 
LboI - Atfhalo relation. 



back from the central AGN can extend up to a few hundred 
kpc and thus can potentially affect the satellite distribu- 
tion. However in higher mass halos satellite AGN will also 
be residing in more massive subhalos and so they will be 
relatively unaffected from the feedback effects of the central 
AGN. We note that in the self regulatory growth paradigm 
the AGN shuts down its own growth by blowing up the gas 
around it (as seen in Fig. 3) but the feedback effects do not 
affect the gas distribution around neighboring black holes. 



3.3 Peak versus Instantaneous Luminosity 

We see a difference in the mean of the CLF for the cen- 
tral AGN between z — 1.0 and z = 3.0. To investi- 
;ate this effect we extracted the light curves of the AGN 
Colberg fc Pi Matted [200j ; lDegraf et al.ll2010l ) and looked 
at their peak luminosities between z — 2.0 and z = 1.0. The 
result is shown in Fig. 7. The peak luminosities (blue stars) 
of the central AGN show a lower scatter with halo mass and 
the best-fit relation is 




Figure 6. Dependence of the conditional luminosity function of 
satellite AGN on central AGN luminosity. The open and the filled 
black circles connected by the dashed and solid lines respectively 
represent the halo mass bin 11.0 < Log(Mhaio) < 12.0. Similarly 
the blue open triangles and the open squares connected by the 
dashed and solid lines respectively represent the halo mass bin 
12.0 < Log(Mhaio) < 13.0. The solid and the dashed lines rep- 
resent the distribution of satellite AGN luminosities for central 
AGN with LboI above (higher) and below (lower) the mean lu- 
minosity (see Fig. 5) respectively. We do not see any significant 
correlation between central AGN luminosity and the number of 
satellites. The error bars at each bin represent the Poisson error 
bars. 
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Figure 7. Bolomctric luminosity of the central AGN as a function 
of halo mass calculated at z = 1.0. The blue points show the 
peak luminosity between z = 2.0 and z = 1.0, and the red open 
circles show the instantaneous luminosity at z = 1.0. The peak 
luminosity tends to correlate more tightly with halo mass than 
instantaneous luminosity. 
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Figure 8. The mean occupation distribution of the AGN as a function of halo mass. The top, middle, and bottom panels correspond 
to redshifts 1.0, 3.0, and 5.0, respectively. The left, middle and right columns correspond to different luminosity cuts. The black points 
represent the mean occupation of all the AGN within R200 an d the red open circles show the contribution from the satellite AGN. The 
mean occupancy of the central AGN can be idealized as a softened step function while the satellite population can be approximated by a 
power law (see Eqs. 4 and 5). The black solid, blue dot-dashed, and the red dashed lines are the best-fit models for the total, central, and 
satellite occupation respectively. The best-fit HOD parameters are shown in Table 2. The error bars reflect the 1<t Poisson error bars. 



where I/Boi is the peak bolometric luminosity. We over- 
plot the instantaneous luminosity at z = 1.0 (red open cir- 
cles). The best-fit slope for the instantaneous luminosity is 
(0.95 ±0.12). As we go to higher redshifts the AGN have lu- 
minosities closer to their peak luminosities and hence we see 
a higher mean value and a tighter correlation (Fig. 5). The 
difference between the peak and the instantaneous luminos- 
ity is most prominent at higher mass halos. We associate this 
suppression with feedback from AGN. In Paper I, we used 
the Mbh — a relation to show that black holes residing in 
higher mass halos enter the feedback dominated phase at an 
earlier time than black holes populating lower mass halos. 
Thus feedback effects will alter the M^aio — £bo1 correlation 
resulting in a wider distribution of host halo masses for a 
given AGN lumin osity at lower red shift. This result is also 
in agreement with lLidz et alj (|200fJ ) who conclude that the 
peak luminosity is more correlated with halo mass than the 
instantaneous luminosity of AGN and hence a better indi- 
cator of clustering. However, we cannot measure the peak 
luminosity of the AGN in practice. 



3.4 The Mean Occupation Function of AGN 

We model the mean occupation function of AGN in 
dark matter halos by decomposing it into a more physi- 
cally illuminating central and satellite contributions (e.g ., 
iKravtsov et~al]|2004l ; IZheng et al.1l2005l ; IZehavi et af]|2005h 

{N (M)> = (AT ccn (M)> + (AU (M)) (3) 

where (N ceu (M)) represents the mean occupation function 
of central AGN and (iV sa t (M)) represents the mean occu- 
pation function of satellites (see discussions in §2.2 for the 
differences in terminology with Paper I) . This formalism im- 
plicitly assumes that the AGN content of halos of a given 
mass is statistically independent of the large-scale environ- 
ments within which the halos reside and thus the mean occu- 
pation (N (M)) depends only on ha lo mass fe.g. jBond et alj 
Il99ll ; iLemson fc Kauffmannlll999l ). (N (M)) is shown as a 
function of halo mass in Fig. 8. The black filled circles show 
the total occupation and the red open circles show the satel- 
lite occupation. The black solid, blue dot-dashed, and the 
red dashed lines are the best-fit models for the total, cen- 
tral, and satellite occupation. The top, middle, and bottom 
panels show redshifts 1.0, 3.0, and 5.0, respectively, while the 
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left, middle, and right columns denote luminosity thresholds 
10 38 ergs/s, 10 40 ergs/s, and 10 42 ergs/s. 

We see that the central AGN occupation number fol- 
lows a distribution close to a softened step function and the 
satellite occupation follows a power la w similar to the galaxy 
or the dark matte r subhalo case (e.g.. iKravtsov et alj|2004c 
IZheng et al1l2005l ). Our HOD model is defined as follows: 



1 + erf 



LogM - LogMmir 



OLogM 

<iVsat) = (M/Mi) a exp(-M cu t/M) 



(4) 
(5) 



In this formalism there are four parameters for modeling the 
HOD: M m i n , defining the halo mass where the occupation of 
central AGN of a given type (in the present case we choose 
AGN in terms of luminosity type) is 0.5; OLogM the char- 
acteristic transition width; Mi, the mass scale at which the 
mean number of satellites above a given luminosity threshold 
equals unity; and a, the power law exponent for the satellite 
occupation. However it has been observed that the mean oc- 
cupation of the satellites drops off faster than a power law 
at lower halo mass (e.g., |Zheng et al.ll2005l : iKravtsov et al.l 
12004 IConrov et aLll200fj| ) and we also see this trend in our 
case. We use the parameterization of iTinker et alj |2005l ) to 
model this drop-off. The parameter M cu t is used to model 
the rolling off the power law. So we have a five parame- 
ter HOD model. The best-fit HOD parameters are shown in 
Table 2. The change in the luminosity threshold does affect 
the mean value of the power law exponent but the values 
are consistent within la. The power law exponent also show 
weak evolution (seen in Paper I) with redshift but is consis- 
tent with no evolution within statistical limits. 

We compute the distribution of AGN with respect to the 
mean P(N\{N)). In Paper I we showed that the distribution 
of satellite number (without any mass or luminosity cut) 
follows a Poisson distribution. We have verified that this is 
also the case for AGN in any luminosity threshold sample. 
We further performed a Kolmogorov-Smirnov (KS) test to 
check whether the satellite distribution follows a Poisson 
distribution in our luminosity selected sample. The mean 
P-value that we obtain by performing the KS test over all 
redshifts and over all mass and luminosity bins is 0.92. This 
shows that the null hypothesis is strongly accepted and the 
distribution of the satellites is close to a Poisson distribution. 



3.5 Luminosity versus Mass 

Our model of central occupation is slightly different than Pa- 
per I. In Paper I we do not impose any mass cut and hence 
the fraction of halos (above the threshold mass) containing 
black holes is always unity (see Eqs. 1 and 2 in Paper I). 
In this paper the softening of the step function arises from 
the luminosity based selection. We note that at lower lumi- 
nosities the duty cycle is extremely close to 1.0 at all halo 
massbins and we converge to the model of Paper I. However 
the drop-off from unity is evident in the 10 42 ergs/s sample 
(right panel of Fig. 8). In Fig. 9 we show the differences in 
the mean occupation function between a mass selected and 
a luminosity selected sample. The black solid lines show the 
mean occupation function of AGN with bolometric luminosi- 
ties greater than 10 40 ergs/s at z — 1.0 (shown in top middle 
panel of Fig. 8). The triangles represent the mean occupa- 
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Figure 9. The mean occupation distribution of the AGN for dif- 
ferent selections. The black solid lines represent the mean occupa- 
tion function of AGN with bolometric luminosities greater than 
10 40 ergs/s at z = 1.0 (shown in the top middle panel of Fig. 8). 
The triangles represent the mean occupation of central AGN and 
the open circles represent the mean occupation of satellite AGN. 
We now use the best-fit re lation between black h ole mass and 
AGN luminosity at z = 1.0 ( Chatterjee ct al. 20Q3) to obtain the 
equivalent mass corresponding to our bolometric luminosity of 
10 40 ergs/s. The red dashed lines represent the mean occupation 
functions for the mass selected sample with triangles represent- 
ing central black holes and open circles representing satellite black 
holes. In all cases the error bars are lcr Poisson error bars. 



tion of central AGN and the open circles represent the mean 
occupation of satellite AGN. We now use the best-fit rela- 
tion between black hole mass and AGN luminosity at z = 1.0 
l|Chatteriee et al.ll2008h to obtain the equivalent mass cor- 
responding to our bolometric luminosity of 10 40 ergs/s. The 
red dashed lines represent the mean occupation functions 
for the mass selected sample with triangles representing cen- 
tral black holes and open circles representing satellite black 
holes. 

We see a clear difference between the two populations at 
lower halo masses. For the central AGN (triangles) the solid 
and the dashed lines converge at a mass scale of 10 11 M©. 
Also the mass-selected occupation function falls-off very 
steeply below this mass scale. This difference arises from 
lower mass black holes residing in lower mass halos with high 
Eddington ratios. Also the correspondence between black 
hole mass and host halo mass is tighter than luminosity and 
hence we see this sharp cut-off in t he occupation funct ion for 
the mass selected sample. Recentlv lGallo et al.l (|2010T ) found 
evidence of such low mass high accretion rate black holes in 
local early type galaxies in the AMUSE- Virgo survey. For 
the satellite population (open circles) the lack of conver- 
gence in the occupation function between the two popula- 
tion is prominent even at higher halo masses. The minimum 
halo mass for hosting satellite black holes (based on mass) is 
much higher than minimum halo mass for hosting satellite 
AGN with equivalent luminosity. We thus see that the HOD 
properties will be significantly different between a mass se- 
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10.41 ±0.06 
10.05 ± 0.08 
11.85 ± 0.05 


0.50 ±0.05 
0.78 ±0.13 
0.59 ±0.06 


0.55 ±0.13 
0.54 ±0.13 
1.00 ±0.32 


11.36 ± 1.48 
11.39 ± 1.41 
12.66 ± 2.12 


11.70 
11.70 
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5.0 


> 10 38 
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> 10 42 
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121 


66 
66 
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10.91 ± 0.04 
10.41 ± 0.06 
11.53 ± 0.09 


0.16 ±0.03 
0.57 ±0.05 
0.52 ±0.18 


0.75 ± 0.18 
0.75 ± 0.18 
1.37 ±0.37 


11.47 ±2.81 
11.47 ±2.81 
12.05 ±3.83 


11.70 
11.70 
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Table 2. AGN HOD parameters for three redshifts corresponding to Eqs. 4 and 5. Columns 3 and 4 show the total number of central 
and satellite black holes in each bin. 



lected sample and a luminosity selected sample. This is be- 
cause a relatively small halo can host a bright AGN if there 
is high density gas. On the other hand, it is difficult to get a 
massive black hole in a small halo, since that would require 
significant amounts of dense gas over a long period of time. 

3.6 Radial Profiles 

Fig. 9 shows the radial distribution of the number density of 
satellite AGN (luminosity greater than 10 40 ergs/s) within 
host halos of masses 11.0 < Log(A/haio) < 13.0. The red cir- 
cles, blue triangles, and the green squares show the profile at 
z = 1.0, z = 3.0, and z = 5.0 respectively. It has been shown 
in Paper I that the radial distribution of black holes follow 
a power law at all redshifts (Eq. 6 of Paper I). Our findings 
are similar to Paper I. The black dashed line is the best- 
fit power law to the average profile over all redshifts (since 
we do not see any significant evolution with redshift). The 
best-fit values for f3 (power law index of AGN: averaged over 
all redshifts) and Log(no) (normalization: averaged over all 
redshifts) are -2.33 ± 0.08 and -0.67 ± 0.05 for the sam- 
ple shown in Fig. 9. The power law index j3 does not show 
any significant dependence on halo mass or AGN luminos- 
ity either. For comparison we also show the profile of dark 
matter and galaxies within AGN host halos. The solid and 
the dotted lines show the average profiles (averaged over all 
redshifts) for dark matter and satellite galaxies respectively. 
The minimum stellar mass that we used for obtaining the 
galaxy profile is 1O 8 M0 (roughly 100 times the stellar mass 
resolution). The profiles do not show any variation if we 
change the threshold stellar mass for selecting galaxies. The 
current choice of stellar mass is an optimization between res- 
olution elements and statistics. The profiles are normalized 
to the mean number density of the corresponding species 
(sate llite AGN/dark matter/ satellite galaxies) within -R200 

fe.g-. lNagai fc Kravtsovll2005l ). 

We also fit NFW profile (|Navarro et all 119971 ) to the 
AGN radial distribution. We fit the profiles for different 
concentration parameters and calculate the corresponding P 
values. At all redshifts our data strongly disfavors the null 
hypothesis and the NFW profile is ruled out at 3a. We note 
that the AGN are more centrally concentrated than dark 
matter and galaxies. The reason that we see an enhanced 
population of AGN at the center of the halo compared to 
galaxies is because of the merging process of black holes. 



When two galaxies merge there exists a time lag between 
merging galaxies and the merging of AGN that reside in 
them. This time lag is due to the time it takes for the satel- 
lite AGN to fall in to the halo center where it can merge with 
the central AGN. After the two galaxies merge the time that 
it takes for the AGN to merge can then be furthe r affec ted by 
the gas content of these galaxies. iLin fc Mohrl (|2007l ) mea- 
sure the radial profile of radio sources in clusters and show 
that it is consistent with an N FW profile with a concentra- 
tion of 25. lMartini et alj l|2007l ) study the radial distribution 
of X-ray selected AGN in clusters and find that AGN with 
X-ray luminosities above 10 42 ergs/s show stronger central 
concentration than cluster host galaxies. However, a bigger 
sample with Lx > 10 41 ergs/s (closer to our sample of AGN) 
does not show any stronger evidence of central concentra- 
tion; different from what we observe in simulation. 

In Fig. 11 we show the relation between stellar mass and 
AGN luminosity at z = 1.0. The stellar mass is computed 
within a spherical region of radius 30 kpc surrounding the 
black hole. We selected different spatial scales to compute 
the stellar mass. The values converge between radii of 25 
kpc to 30 kpc and hence we chose 30 kpc to be the relevant 
radius for computing the stellar mass content in AGN host 
galaxies. The black points represent central AGN and the 
red open circles show the satellite distribution. In the case 
of central AGN the stellar mass for more luminous AGN 
is generally higher. AGN will have higher luminosity when 
the accretion rate is high, which is more likely when the 
gas density is high. Also when the gas density is high the 
cooling rate will be higher and hence there will be more 
stars. So in general higher density will correspond to high 
accretion rates of AGN and higher stellar mass. For satellite 
AGN this trend is very weak. But even for satellite AGN 
if we consider a luminosity selected sample we will system- 
atically choose AGN with higher stellar mass. Also in Fig. 
3 we see that AGN with higher stellar mass content tend 
to reside in the central regions of the halos. This can po- 
tentially increase the higher concentration of AGN seen in 
Fig. 10 for a luminosity selected sample. To further check 
this we selected AGN based on their stellar mass content 
and computed their average radial profile. We find a steeper 
mean slope (f3 ~ —2.7) for the stellar mass selected sam- 
ple compared to the luminosity selected sample. This again 
is an indication of enhanced central concentration of AGN 
with higher stellar mass content in their host galaxies. 
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Figure 10. The radial profile of satellite AGN. The open circles, 
triangles, and squares show the profiles at z = 1.0, z = 3.0, and 
z = 5.0 respectively. The profiles are obtained by averaging over 
AGN with LboI > 10 40 ergs/s and host halo masses between 
10 11 — 1O 13 M0. The radial profiles do not exhibit any significant 
evolution with redshift. The dashed line shows the average power- 
law fit over all redshifts. The best-fit value for the slope is — 2.33± 
0.08. For comparison we also show the profile of dark matter 
and galaxies within AGN host halos. The solid and the dotted 
lines show the average profiles (averaged over all redshifts) for 
dark matter and satellite galaxies respectively. The error bars are 
representative of the ltr error bars. 



4 SUMMARY AND DISCUSSIONS 

In this paper, we study the relation between AGN and dark 
matter halos using a cosmological SPH simulation that in- 
corporates black hole and galaxy formation physics. Specifi- 
cally we have investigated the effect of environmental factors 
(e.g., feedback, local gas density, host galaxy mass) on the 
accretion rates of AGN within dark matter halos and how 
it affects the occupation distribution of AGN. We examine 
a number of simplifying assumptions in the HOD modeling 
(e.g., dependence of AGN luminosity on halo mass, depen- 
dence of AGN accretion on black hole mass) of AGN and 
provide the necessary tools to model AGN clustering data. 
We have characterized the HOD of faint AGN (the sample 
probed by our simulation box) which can be generalized to 
incorporate the bright end of the luminosity function. 

We compute the conditional luminosity functions of 
AGN and separate the contributions from the central and 
satellite AGN. Our key findings are as follows. (1) The cen- 
tral AGN luminosities follow a log-normal distribution simi- 
lar to the assumption in semi-analytic studies. (2) The mean 
of the CLF for a given halo mass shows a strong dependence 
on redshift with higher luminosities at higher redshifts. (3) 
The scatter in central AGN luminosity is large at low red- 
shift, but decreases with increasing redshift. This implies 
that the dependence of AGN clustering on luminosity is 
weak at low redshift, but it can become stronger at high 
redshift. We analyze the light curves of individual AGN and 
show that there exists a tighter correlation between halo 



mass and peak luminosity rather than instantaneous lumi- 
nosity. We present the joint distribution of satellite occupa- 
tion as a function of halo mass and central AGN luminos- 
ity. We do not see any significant correlation between the 
satellite number and the luminosity of the central AGN. We 
also show that the mean occupation function of the central 
AGN resembles a softened step function while the satellite 
population follows a power law with an exponential roll-off 
at lower mass, similar to what has been observed with the 
galaxy HOD. We show that low mass black holes with high 
Eddington ratios residing in low mass halos makes the lumi- 
nosity based HOD significantly different from the black hole 
mass based HOD. 

We will now compare our simulation results with semi- 
analytic models that have been widely used for cluster- 
ing an alysis (e.g.. [Martini fc W einberg 2001; Shanka r et al.l 
l2010ah. We compare our results with the model described 
in lShankar et ail f2010a) who computed the mass function 
of black holes at z — 3.0 using the halo mass function and 
a redshift dependent relation between halo mass and black 
hole mass. They parameterize the redshift dependence in 
the normalization of the black hole mass-halo mass relation 
while the slope is kept fixed at 1.52. In our simulation we ob- 
serve the average slope to be (1.04 ±0.06) at z = 3.0. There 
is a large scatter from the mean slope and the la scatter is 

0. 44. We note that this relation in our simulation is depen- 
dent on the AGN model and the ratio of the threshold halo 
mass that can host b l ack hol es to the mass of the seed black 
holes. IShankar et al.l l|2010ah assumes the black holes to be 
accreting with a constant Eddington fraction. Although the 
mean \ in our simulation is roughly independent of halo 
mass as seen in Figure 1, they do s how an evolution with 
redshift. This has been also noted in IShankar et all (|2010aT ) 
where they suggest that the assumption of constant Edding- 
ton fraction might break down at the faint end. The other 
important parameter is the scatter in the LboI — A^haio rela- 
tion. We find a log-normal distribution as assumed in semi- 
analytic studies. However the scatter strongly depends on 
redshift with increasing scatter at low redshifts (shown in 
Figure 5). 

Observationally the most relevant physical quant i ty de - 
scri bing an AGN is l u minos ity. Recentlv lMivaii et al.l (|2010h 
and IStarikova et al.l (|2010h have used X-ray selected AGN 
to constrain the HOD empirically. We propose to compare 
our HOD model with observational samples in a follow-up 
paper. However an alternative approach would be select- 
ing black holes based on their mass and measure clustering 
statistics based on black hole mass. Although we do not have 
reliable observational measurements of black hole masses, we 
have also provided the framework for predicting clustering 
properties based on the black hole mass function in Paper 

1. These two complimentary approaches would impose even 
tighter constraints on theoretical models of AGN growth and 
feedback. Our analysis provides the first step toward com- 
paring semi-analytic and simulation results on AGN clus- 
tering and assessing some of the simplifying assumptions in 
the present interpretations of AGN clustering observations 
and constraining physical parameters. Because of the small 
simulation box, our current study is limited to low luminos- 
ity AGN, and the results on satellite AGN suffer from small 
number statistics. A larger simulation box with larger sta- 
tistical samples spanning both the faint and the bright end 
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Figure 11. Relation between stellar mass and AGN luminosity at 
z = 1.0. The stellar mass is computed within a spherical region 
of radius 30 kpc surrounding the black hole. The black points 
represent central AGN and the red open circles show the satellite 
distribution. In the case of central AGN the stellar mass for more 
luminous AGN is generally higher. For satellite AGN this trend is 
very weak. But even for satellite AGN if we consider a luminosity 
selected sample we will systematically choose AGN with higher 
stellar mass. 



of the luminosity function can provide tighter constraints 
on AGN HOD. Also, the mass of the seed black hole is de- 
pendent on the resolution of the simulation, imposing an 
artificial mass cut in the simulation. Future work on high 
resolution simulation an d more accurate mode ling of accre- 
tion and feedback (e.g., iBooth fc Schavell2009l ) will also be 
needed to understand the full implication of gas physics and 
black hole accretion on clustering studies. The HOD formal- 
ism has been successfully incorporated in galaxy evolution. 
We hope that our work will have the same impact for study- 
ing the co-evolution of AGN and their hosts and shed light 
on the AGN-galaxy connection. 
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